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LONG-TERM  GOAL 

The  long-term  objective  is  to  develop  noninvasive  instrumentation  and  techniques  for  in  situ 
measurements  of  the  interfacial  properties  of  the  ocean,  including  the  intrinsic  viscoelastic 
properties.  These  intrinsic  properties,  as  interpreted  through  the  laws  of  mechanics,  will  be 
related  to  the  apparent  properties  of  the  interface  as  deduced  from  the  usual  phenomenological 
approach  of  short-wave  dampening.  Ultimately,  the  intrinsic  interfacial  properties  are  required 
for  computational  models  of  the  ocean  surface,  derived  from  first  principals,  to  predict  the  shape, 
motion  and  other  factors  that  influence  radar  scatter  from  a  surfactant-influenced  interface  such 
as  the  ocean  surface. 

The  rapid  increase  in  computing  power  in  the  past  two  decades  has  made  it  possible  to 
numerically  solve  the  exact  equations  of  fluid  motion  over  a  moderate  range  of  scales  in  time  and 
space.  And  although  the  equations  of  motion  for  the  flow  in  the  air  and  in  the  sea,  namely  the 
Navier-Stokes  equations,  have  been  known  for  over  a  century  and  a  half,  the  boundary 
conditions  for  these  equations  are  generally  unknown  at  a  surfactant-influenced  gas/hquid 
interface  because  the  boundary  conditions  depend  on  the  intrinsic  (cf  apparent)  viscoelastic 
properties  of  the  interface,  measurements  of  which  were  not  feasible  until  recently 

SCIENTIFIC  OBJECTIVES 

The  scientific  objective  of  the  present  project  is  to  determine  the  intrinsic  viscoelastic  properties 
that  have  a  significant  effect  on  the  movement,  shape,  and  other  factors  that  influence  radar 
scatter  from  a  surfactant-influenced  gas/hquid  interface  such  as  the  ocean  surface. 

InitiaUy,  through  the  use  of  the  most  general  stress-strain  relation  for  a  Newtonian  interface, 
namely  the  Boussinesq-Scriven  constitutive  relation,  we  wish  to  directly  determine  the 
viscoelastic  properties  of  the  interface.  This  must  first  be  performed  in  the  laboratory  for  weU- 
behaved  surfactant  systems  and  well-controUed  subsurface  flows  and  waves,  prior  to  the 
application  of  the  technique  to  in  situ  measurements.  Ultimately,  non-Newtonian  constitutive 
relations  wiU  be  utilized  for  the  evaluation  of  the  viscoelastic  properties  of  natural  polymeric 
films  on  the  ocean  surface. 

There  are  two  coefficients  in  the  Boussinesq-Scriven  that  describe  the  dissipative  or  the  viscous 
part  of  the  stress  at  the  interface.  These  are  the  interfacial  shear  viscosity  (|i‘^)  and  interfacial 

dilatational  viscosity  (K'^).  In  general,  these  coefficients  depend  on  the  thermodynamic  state  of 
the  surfactant-influenced  interface,  i.e.  the  surface  concentration  of  the  surfactant.  In  addition  to 
the  viscous  effect  at  surface,  the  elasticity  of  the  surface,  which  results  from  gradients  in  the 


1 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1.  REPORT  DATE 

30  SEP  1997  ‘REPORT  TYPE 

3.  DATES  COVERED 

00-00-1997  to  00-00-1997 

4.  TITLE  AND  SUBTITLE 

Interaction  of  Surfactants  with  Shear  Flows  and  Surface  Waves 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Rensselaer  Polytechnic  Institute, Department  of  Mechanical  Engineering, 
Aeronautical  Engineering  and  Mechanics, Troy, NY, 12180-3590 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  aS 

unclassified  unclassified  unclassified  Report  (SAR) 

7 

Standard  Form  298  (Rev.  8-98} 

Prescribed  by  ANSI  Std  Z39-18 


Office  of  Naval  Research 


Space  and  Remote  Sensing 
1997  Annual  Report 

thermodynamic  surface  tension  must  be  quantified.  The  fluid  motion  adjacent  to  the  interface  is 
coupled  to  the  interface  due  to  the  viscosity  in  the  bulk  of  the  fluid.  Thus,  in  order  to  quantify  the 
viscoelastic  effect  at  the  interface,  the  bulk  fluid  motion  adjacent  to  the  interface  as  weU  as  the 
composition  and  motion  of  the  interface  itself  must  be  measured.  Reproducible  and  consistent 
measurements  of  the  intrinsic  viscosities,  especially  are  yet  to  appear  in  the  hterature  for  any 
surfactant  system.  Pubhshed  measurements  of  for  a  given  surfactant  vary  by  as  much  as  a 
factor  of  10  ,  depending  on  the  method  used  to  measure  k  . 

APPROACH 

Well-controlled  flow  fields  are  being  utihzed  in  the  laboratory  in  order  to  study  the  viscoelastic 
properties  of  a  water  surface  covered  by  various  monolayers.  The  unsteady  dilation  and 
compression  of  the  ocean  surface  resulting  from  atmospheric  turbulence  (wind),  shear  in  the 
hquid  side  and  surface  or  internal  waves  is  simulated  in  the  laboratory  by  the  flow  field  due  to  the 
rise  and  consequent  interaction  of  a  laminar  vortex  pair  (see  Figure  1).  The  vortex  pair  stretches 
the  surface  without  intrusive  barriers,  such  as  in  Langmuir  trough-type  of  devices.  Also,  we  have 
recently  developed  other  flow  facilities  including  a  c-g  wave  tank  system  as  weU  as  a  steady-state 
flow  loop  with  a  film  under  compression  by  a  uniform  subsurface  flow  (“Reynolds  ridge  device”) 
designed  primarily  as  a  benchmark  test  for  the  new  areal  SHG  system.  The  boundary  layer 
beneath  the  surface  is  examined  in  each  system.  The  c-g  waves  are  used  as  an  extension  of  the 
vortex  pair  study  in  order  to  stretch  and  dilate  the  surface,  thus  allowing  the  measurement  of  the 
viscoelastic  properties  of  the  surface.  The  c-g  waves  are  also  used  to  study  the  relationship 
between  the  intrinsic  viscoelastic  properties  and  wave  dampening  characteristics,  which  affect 
the  radar  back  scatter  from  the  ocean  surface. 

The  velocity  measurements  in  the  bulk  of  the  fluid  must  include  the  free  surface  boundary  layer, 
which  for  high  Reynolds  number  flows  can  be  very  thin  (of  the  order  of  a  milhmeter  or  less). 
Thus,  the  noninvasive  laser-based  method  of  digital  particle  image  velocimetry  (DPIV)  is 
undergoing  further  refinement  in  our  laboratory  to  improve  the  spatial  and  temporal  resolution  of 
the  velocity  field  measurements.  In  order  to  comphment  the  surface  viscosity  measurement 
techniques  that  we  have  developed,  a  deep-channel  surface  viscometer  is  also  being  tested.  The 
deep-channel  viscometer  provides  measurements  of  the  interfacial  shear  viscosity  (|i‘^)  which  can 
be  subtracted  from  the  measured  value  of  the  sum  of  the  surface  viscosities  to  provide  K'^. 

WORK  COMPLETED 

Many  tasks  have  been  accomplished  in  the  first  year  of  the  project.  A  new  laboratory  has  been 
estabhshed.  The  new  laboratory  is  located  in  a  building  with  an  underground  connection  to  Prof. 
Korenow ski's  laboratory.  This  has  permitted  safe  transportation  of  dehcate  equipment,  including 
the  high  frequency  (100  Hz)  Infinity  Laser  and  support  systems,  to  and  from  each  laboratory 
without  exposure  to  the  elements.  Our  new  laboratory  has  been  equipped  with  a  high  capacity 
water  distillation  system  and  air  filtration  facihties.  A  large  clean-water  tank  with  a  600  liter 
capacity  for  vortex  measurements  has  been  put  into  operation.  The  deep-channel  viscometer  was 
completed  and  put  into  operation. 

A  Reynolds  ridge  device  (see  Figure  2)  was  constructed  of  inert  materials  (glass,  acrylic  plastic, 
and  stainless  steel)  and  is  being  used  initially  as  a  tool  to  tests  and  calibrate  the  areal  SHG  system 
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in  flowing  systems.  Since  the  strain  and  the  resulting  stress  in  the  boundary  layer  beneath  the  film 
can  be  measured  directly  with  laser  velocimetry,  the  surfactant  concentration  gradient  on  the 
surface  can  be  determined  and  used  for  evaluating  the  concentration  profile  measured  by  the 
SHG  camera. 

RESULTS 

The  instrumentation  for  noninvasive  probing  of  surfactant  films  in  the  laboratory  (via  SHG  and 
DPIV)  was  developed  and  the  results  have  been  pubhshed  (Hirsa,  et  al.  1997a).  Our 
measurements  provide  data  on  the  evolution  of  the  surfactant  concentration  field  during  the 
interaction  of  a  transient  subsurface  flow.  The  analysis  presented  in  that  paper  showed,  for  the 
first  time,  that  quantitative  information  on  the  flow  in  the  bulk  fluid  can  be  directly  obtained  from 
nonhnear  optical  probing  of  the  surface  composition.  Also,  a  novel  technique  for  the  direct 

determination  of  the  sum  of  the  surface  dilatational  and  shear  viscosities  (K'^  +  |i‘^)  using  SHG  and 
DPIV  data  was  developed  for  insoluble  surfactants  and  the  results  were  published  (Hirsa,  et  al. 
1997b).  The  extension  of  this  technique  to  soluble  systems,  and  ultimately  to  natural  ocean  films 
is  the  central  thrust  of  the  present  effort  and  is  expected  to  become  feasible  once  the  areal  SHG 
system  is  fully  tested.  Toward  that  end,  a  Reynolds  ridge  device  was  built  and  tested  and  the 
results  foUow. 

A  new  flow  loop  was  specially  designed  and  constructed  to  study  surfactant  films  under  steady- 
state  compression  and  to  provide  a  benchmark  test  of  the  areal  SHG  system  in  flowing  systems 
(see  Figure  2).  An  existing  Laser-Doppler  Velocimeter  (LDV)  system  has  been  modified  to 
provide  scanning  velocity  measurements  with  high  spatial  resolution  (50  microns)  simultaneously 
with  the  areal  SHG  measurements.  The  LDV  measurements  (i.e.  pointwise  measurements  made 
relatively  quickly  along  the  surface-normal  direction)  provide  the  instantaneous  shear  stress  on 
the  film  which  must  be  balanced  by  the  gradient  in  the  surface  tension  and  hence  surfactant 
concentration.  Prehminary  LDV  measurements  of  the  shear  stress  beneath  a  hemicyanine 
monolayer  film,  made  in  the  new  flow  loop,  are  presented  in  Figure  3,  along  with  comparisons  to 
the  laminar  boundary  layer  theory  on  a  sohd  wall.  Figure  3  shows  that  although  the  stress 
distribution  beneath  the  surface  film  is  hnear  when  plotted  with  the  appropriate  scaling  in 
agreement  with  the  theory  for  a  solid  wall;  interestingly,  the  stress  is  greater  than  the  stress 
beneath  a  sohd  wall.  The  SHG  data  that  was  acquired  simultaneously  along  with  the  above  LDV 
measurements  are  expected  to  be  analyzed  in  the  next  few  days  and  will  be  presented  at  the  APS 
meeting  in  San  Francisco  (Korenowski,  Hirsa,  Van  Wagenen  &  Jin  1997). 

A  novel  technique  has  been  developed  and  successfully  tested  which  provides  the  instantaneous 
shape  of  the  air/water  interface,  thus  permitting  the  DPIV  measurements  to  be  performed  in 
boundary -fitted  coordinates  (i.e.,  measurements  on  grids  fitted  to  the  deforming  air/water 
interface  instead  of  orthonormal  coordinates  used  in  conventional  DPIV).  In  the  boundary-fitted 
DPIV  technique,  the  particle  images  are  first  utilized  to  map  the  location  of  the  surface  through 
an  robust,  automated  cross-correlation  method.  This  is  possible  because  of  the  total  internal 
reflection  that  occurs  when  viewing  submerged  objects  from  below  the  surface.  Surface 
elevations  of  up  to  2  mm  at  any  x-position  can  thus  be  determined  to  an  accuracy  of  better  than 
10  microns.  This  technique  provides  a  two-dimensional  mapping  of  the  surface.  The  boundary- 
fitted  measurements  are  particularly  important  because  of  the  large  gradients  that  exist  in  the 
velocity  within  the  free-surface  boundary  layer.  The  results  of  this  work  will  be  presented  at  the 
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APS  meeting  (Mortzheim,  Gayton  &  Hirsa  1997).  Representative  measurements  with  the 
boundary  fitted  DPIV  are  presented  in  Figure  4.  The  variation  of  the  surface  location  p  (averaged 
over  the  4  mm  field  of  view)  as  a  function  of  time  is  shown  in  Figure  4(a).  Figure  4  shows  the 
surface  upwelhng  due  to  the  arrival  of  the  vortex  pair.  The  surface  deformations  associated  with 
the  secondary  vortex  structure  is  also  evident  in  Figure  4.  The  surface  velocity  us  is  presented  in 
Figure  4(b)  and  its  surface-normal  gradient,  i.e.  the  surface  viscosity,  is  shown  in  Figure  4(c). 
Figure  4  shows  that  the  vortex  pair  flow  field  dilates  the  surface  until  the  surfactant  is  completely 
cleared  away  from  the  center  for  a  period,  this  process  continues  until  the  relatively  larger 
surface  tension  of  the  clean  region  pulls  back  the  film  toward  the  center. 

IMPACT/APPLICATION 

We  have  shown  for  the  first  time  that  the  intrinsic  viscoelastic  properties  of  surf  actant- influenced 
surfaces  may  be  measured  directly  through  optical  techniques.  These  intrinsic  properties  are 
required  in  numerical  simulations  of  fluid  motion  on  each  side  of  the  air/sea  interface. 

TRANSITIONS 

We  expect  to  develop  the  instrumentation  and  techniques  to  perform  boundary-fitted  DPIV 
measurements  in  the  free  surface  boundary  layer  beneath  ambient  waves  on  the  ocean  within  the 
next  two  years;  this  will  compliment  the  nonhnear  optical  measurements  of  the  surface. 

RELATED  PROJECTS 

As  described  throughout  this  report,  the  present  research  project  is  closely  linked  to  the  research 
in  nonlinear  optical  probing  of  the  ocean  surface,  conducted  by  Prof.  Korenowski  in  the 
Chemistry  Department  at  RPI.  Our  synergistic  approach  functions  in  such  a  way  that  their  efforts 
will  continue  to  provide  us  with  the  complementary  tools  needed  for  measurements  of  the 
viscoelastic  interfacial  properties  in  the  laboratory  and  ultimately  in  situ,  and  our  flow  facihties 
and  flow  measurements  provide  the  benchmark  for  the  development  and  testing  of  their  surface 
probes  prior  to  their  implementation  for  in  situ  studies. 
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Figure  3.  Shear  Stress  Beneath  Surfaee  Film,  'to  as  Measured  by  the  LDV. 
Note:  The  sohd  wall  solution  is  due  to  Blausius. 
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Figure  4.  Representative  measurements  with  the  boundary-fitted  DPIV  technique. 

It  is  showing  the  [spatially-averaged]  surface  location,  h,  as  a  function  of  time  (a); 
interfacial  velocity,  Us  (b);  and  [surface-normal]  shearing  strain  rate  [=  surface  vorticity] 
du/dy  (c).  These  measurements  were  made  at  x  =  4  cm  and  n  =  3.2  dynes/cm. 
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